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Abstract 

In this paper, forced convection heat transfer of nanofluids 
over a horizontal flat plate embedded in a porous medium 
saturated with a nanofluid is numerically analyzed. In the 
boundary layer, heat can be generated or absorbed. It is 
assumed that the nanoparticles are uniformly dispersed in 
the base fluid. A similarity approach is used to reduce the 
governing partial differential equation to an ordinary 
differential equation. The resulting ordinary differential 
equation is numerically solved for a type of porous medium, 
sand, and three types of nanoparticles, namely, alumina 
(AI2O3), copper (Cm), and titanium dioxide (T1O2). The effect 
of heat generation/absorption as well as volume fraction of 
nanoparticles on the heat transfer enhancement of 
nanofluids is theoretically investigated. 
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Introduction 

Nowadays, nanotechnology provides a new 
opportunity to produce nanoparticles with dimensions 
below 50 nm (Wang and Mujumdar 2008). Dispersing 
nano-size solid particles in a conventional heat 
transfer fluid would increase the thermal conductivity 
of the fluid. Nanofluid is a fluid which contains 
nanometer-sized solid particles. Successful 
employment of nanofluids reduces the size of heat 
exchanger systems and also enhances the performance 
of these systems (Wang and Mujumdar 2008). 
Nanofluids have many industrial applications such as 
nuclear reactors, transportation, electronics as well as 
biomedicine and food (Kakac and Pramaumjaroenkij 
2009). 



Forced convection heat transfer in porous media also 
has many applications such as geothermal energy 
engineering, groundwater pollution transport, 
chemical reactors engineering, nuclear waste disposal, 
insulation of buildings and pipes (Zhao and Song 
2001). 

Forced convection heat transfer over a horizontal plate 
embedded in a porous medium has been studied in 
many literatures. Beckermann and Viskanta (1986) 
studied forced convection boundary layer flow and 
heat transfer along a flat plate embedded in a porous 
medium. Postelnicu et al. (2001) discussed the effect of 
variable viscosity on the forced convection flow and 
heat transfer past a horizontal surface placed in a 
saturated porous medium in the presence of internal 
heat generation. 

Noghrehabadi et al. examined the flow and heat 
transfer of nanofluids over an isothermal flat plate 
with partial slip boundary condition (Noghrehabadi et 
al. 2012c), suction or injection (Noghrehabadi et al. 
2012a) and magnetic field effects (Noghrehabadi et al. 
2012b). Maghrebi et al. (2012) analyzed forced 
convection heat transfer of nanofluids in a porous 
channel. Rosea et al. (2012) investigated the non-darcy 
mixed convection of nanofluid over a horizontal 
surface placed in a porous medium. 

Forced convection heat transfer of nanofluids over a 
horizontal flat plate embedded in a saturated porous 
medium in the presence of heat generation/ absorption 
has not been analyzed yet. The present study aims to 
analyze the effect of heat generation or absorption on 
the forced convection heat transfer over a horizontal 
plate embedded in a porous medium saturated with a 
nanofluid. The governing partial differential equation 
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is transformed into an ordinary differential equation 
by aid of similarity variables. The governing equation, 
in its non-dimensional form, is a function of heat 
generation parameter and the ratio of thermal 
diffusivity of nanofluid and the base fluid. 

Mathematical Formulation 

Consider the steady two dimensional (x, y) forced 
convection heat transfer of a nanofluid over an 
isothermal horizontal surface which is embedded in a 
saturated porous medium. The schematic view of 
physical and coordinate system is illustrated in FIG. 1. 
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T, 
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Flat plate with constant temperature J",, 

FIG. 1 PHYSICAL MODEL AND COORDINATE SYSTEM 

The nanofluid flows with a constant velocity of U 
through the porous medium because of an external 
pressure gradient. The temperature of the plate is held 
at constant value of Tw. The nanofluid is hemogenius 
and the nanoparticles are uniformly dispersed in the 
base fluid; therefore, the concentration gradient in the 
boundary layer is neglected. 

By assuming an isotropic Darcy porous medium and 
neglecting the wall effects, the governing boundary 
layer equations in the Cartesian coordinate system can 
be written as (Nield and Bejan 2006): 
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the boundary conditions are as follows: 

y = 0: v = T = T W 

y — >• oo : u=U T = T x 



(T-TJ 



(1) 



(2) 



(3) 



(4) 



(5) 
(6) 



where p is the pressure, u is the Darcy velocity and T is 
the temperature of the nanofluid. The effective 
thermal diffusivity of the nanofluid and porous 
medium is denoted by <Xm,nf. Here, K is the permeability 
of the porous medium, \i n ; is the dynamic viscosity of 
the nanofluid. Qo is the heat generation/ absorption 
rate inside the boundary layer. The subscript °° 
denotes the large values outside the boundary layer. 
The thermo physical properties of nanofluid can be 
evaluated as follow: 
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Here, e is the porosity of the porous medium, km,nf is 
the effective thermal conductivity of nanofluid and 
porous medium, km,f is the effective thermal 
conductivity of the base fluid and porous medium, 
(pCp)nf is the heat capacity of the nanofluid, (p is the 
nanoparticle volume fraction, p/ and pp are the 
densities of the base fluid and the nanoparticles, 
respectively, km is the thermal conductivity of porous 
medium, ki is the thermal conductivity of the fluid 
which flows through porous (nanofluid or pure fluid) 
and knf is the thermal conductivity of the nanofluid. kf 
and kp are thermal conductivities of the fluid and the 
nanoparticles, respectively. Eq. (7-c) is based on the 
work of Alizad et al. (2012) and Eq. (7-d) is based on 
the work of Bachok et al. (2011) and Khanafer and 
Vafai (2011). In order to evaluate km,nf, first of all, knf 
should be evaluated from Eq. (7-d) and then the value 
of knf should be used as ki in (7-c). keg is km,f when h is 
the thermal conductivity of pure fluid. 

Eqs. (2) and (3) are simplified using cross- 
differentiation, and the continuity equation is also 
satisfied by introducing a stream function, (ip): 



U = 



dy/ 



v = 



dy/ 

dx 



Solution of Eqs. (2) and (3) leads to: 
u=U v=0 



(8) 



(9) 
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where K , and U is the free stream velocity. In 

order to reduce the governing equation, Eq (4), subject 
to the boundary conditions, Eqs. (5) and (6), to an 
ordinary differential equation, the similarity variables 
of r) and are introduced as, 
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local Peclet number is defined as: 
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Employing the similarity variables, Eq. (10) and (11), 
on the Eq. (4) reduces the governing partial differential 
equation to the following ordinary differential 
equation (See Appendix A): 

a D 0" + ^7]0' + Ae = O (13) 
subject to 

0(0) =1 (14) 
0(oo) = (15) 
where prime denotes differentiation with respect to rj. 
The non-dimensional parameters of old and A are 
defined as follows: 
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Ao is the heat generation/absorption constant, and A is 
the heat generation/absorption parameter, whose 
positive values represent heat generation and 
negative values represent heat absorption. 

In this study, the quantity of practical interest, which 
shows the rate of heat transfer from the plate, is the 
local Nusselt number (Nux): 



k m ,f (T w T m ) 



where 

tfw ~ K m,nf\ )y=Q 

the reduced Nusselt number is written as follow: 
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Results and Discussion 

In order to analyze the effect of nanoparticles on the 
heat transfer, three nanoparticles types, namely, Cu, 
AI2O3, and T1O2 are chosen. The porous medium is 
sand, and the water is the base fluid. The required 
thermophysical and mechanical properties of the 
nanoparticles, porous medium and base fluid are 
shown in TABLES 1 and 2. 

TABLE 1 PROPERTIES OF FLUID AND THREE NANOPARTICLES 
(BACHOK ETAL.2011) 



Physical 
properties 


Fluid 
phase 
(Water) 


Cu 


AI2O3 


TiOz 


C P (J/kgk) 


4179 


385 


765 


686.2 


p(kg/m 3 ) 


997.1 


8933 


3970 


4250 


k (W/m.k) 


0.613 


400 


40 


8.9538 



TABLE 2 PROPERTIES OF SAND POROUS MEDIA (NIELD AND 
BEJAN 2006) 



Physical properties 



Sand 



0.45 
3 



The governing ordinary differential equation, Eq. (13), 
subject to the boundary conditions, Eqs. (14) and (15), 
is solved using Range-Kutta method with shooting 
technique. The asymptotic value of infinity is assumed 
as r/~=8. 

Assuming A=0 (no heat generation/absorption) and 
(p=0 (no nanoparticle dispersed in the base fluid), the 
present work reduces to the work of Bejan (1984). 
Therefore, in order to test the accuracy of the present 
solution, the value of d'(0), obtained by Bejan (1984), is 
compared with d'(0) obtained in the present study in 
the TABLE 3. Moreover in FIG. 2, the temperature 
profile, obtained by Bejan (1984), is compared with 6(r\) 
in the present work when Ao=0 and cj)=0. These 
comparisons show exellent agreement between the 
present study and the work of Bejan (1984). 

The results are obtained for different values of 
nanoparticle concentrations as well as heat generation 
parameter A. 

FIG. 3 shows the effect of heat generation/ absorption 
parameter on the temperature profiles. This figure 
shows that increase of r\ would decrease the local 
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temperature. The temperature on the surface is one 
and decreases to reach its asymptotic value of zero far 
from the plate. It is clear that heat generation gives an 
increase in the temperature of the nanofluid whereas 
heat absorption decreases the temperature of the 
nanofluid. 

Variation of reduced Nusselt number with volume 
fraction of nanoparticles for three types of 
nanoparticles is shown in FIG. 4. In this figure the 
porous medium is sand and the Ao=400000. FIG. 4 
shows that the reduced Nusselt number would 
increase with an increase of nanoparticles volume 
fraction. It should be noted that for a constant value of 
nanoparticle volume fraction, the copper nanoparticles 
provide the highest value of reduced Nusselt number, 
and after that the aluminium oxide nanoparticles 
provide high values of reduced Nusselt number. It is 
worth noticing that the thermal conductivity and heat 
capacity of copper nanoparticles also are higher than 
that of titanium dioxide and aluminum oxide 
nanoparticles. 



TABLE 3 COMPARISON BETWEEN RESULTS OF BEJAN (BEJAN 
1984) AND PRESENT WORK FOR ( X o=0 AND <5=0) 



0'(O) 


Bejan(bej'an 1984) 


Present work 


0.564 


0.564 



3 




























2.5 










* 


Present work 
- Bejan 




2 
















^15 
















1 
















05 




















0.2 


0.4 


e 


0.6 


0.8 1 



FIG. 2 COMPARISON BETWEEN TEMPERATURE PROFILE OF 
PRESENT STUDY AND THAT OF BEJAN (1984) 




FIG. 3 EFFECTS OF HEAT GENERATION PARAMETER ON 
TEMPERATURE PROFILES FOR SAND POROUS MEDIA AND 
AL2O3 NANOPARTICLES WHEN O=0.2 




FIG. 4 VARIATION OF REDUCED NUSSELT NUMBER WITH <P 
FOR SAND AND THREE TYPES OF NANOPARTICLES WHEN 

\ o=400000 

Conclusion 

Effect of heat generation or absorption as well as types 
and concentration of nanofluids on the forced 
convection heat transfer over a horizontal flat plate 
embedded in a porous medium is investigated 
numerically. The results can be summarized as follows: 

1. Heat generation increases temperature of nanofluid 
whereas heat absorption decreases the temperature. 
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2. For all studied cases, an increase of nanoparticles 
volume fraction would significantly increase the 
reduced Nusselt number. 
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Nomenclature 

specific heat at constant pressure C,, 

permeability K 

thermal conductivity k 

local Nusselt number Nu x 

local Peclet number Pe x 

fluid temperature T 

ambient temperature T„ 

velocity components u,v 

Cartesian coordinates x,y 

Greek symbols 

thermal diffusivity a 

porosity £ 

similarity variable 1 

dimensionless temperature 6 

heat generation/absorption parameter X 

dynamic viscosity M 

density p 

nanoparticle volume fraction <P 

Subscripts 

ambient condition 00 

fluid / 

porous media m 
effective value for porous medium and mj 
base fluid 

effective value for nanofluid and porous m , n f 
medium 

nanofluid "/ 

particle P 



Superscript 

differentiation with respect to r] 
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Appendix A: 

Details of transforming the PDE (Energy equation) to 
ODE: 



Differentiation r\ with respect to x and y: 

Si] _ 1 ( U y/2 T - 3 /2 _ Pel' 1 

dx 2 a„ 



(Al) 



dy x 

m,f J 

and differentiation T with respect to x and y: 
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a m,f a m,f 
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Energy equation follows as: 

d 2 T , Go 



rr rr • 



u 1- v — = a 



dx dy m ' nf dy 2 (pc) nf 
Substituting differentiations into energy equation: 
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